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§ ■ ABSTRACT 

MERLIN measurements of 1.6-GHz OH masers associated with Orion-BN/KL are 
\^ , presented, and the data are compared with data on other masers, molecular lines, 

■ compact radio continuum sources and infrared sources in the region. OH masers are 

detected over an area 30 arcsec in diameter, with the majority lying along an approxi- 
mately E-W structure that extends for ~18 arcsec, encompassing the infrared sources 

t> ' IRc2, IRc6 and IRc7. Radial velocities range from -13 to +42 km s _1 . The system of 

OH masers shows a velocity gradient together with non-circular motions. The kine- 

■ matics are modelled in terms of an expanding and rotating disc or torus. The rotation 

axis is found to be in the same direction as the molecular outflow. There is an inner 
cavity of radius ~1300 au with no OH masers. The inner cavity, like the H2O 'shell' 
masers and SiO masers, is centred on radio source I. Some of the OH masers occur 

\ in velocity-coherent strings or arcs that are longer than 5 arcsec (2250 au). One such 

"--»« . feature, Stream A, is a linear structure at position angle ~45°, lying between IRc2 

1-^ ' and BN. We suggest that these masers trace shock fronts, and have appeared, like a 

vapour trail, 200 yr after the passage of the runaway star BN. The radio proper mo- 
Q ■ tions of BN, source I and source n project back to a region near the base of Stream A 

\ that is largely devoid of OH masers. The 1612-MHz masers are kinematically distinct 

c/3 , from the other OH masers. They are also more widely distributed and appear to be 

associated with the outflow as traced by H2O masers and by the 2.12-/j,m emission 
from shocked H2. The magnetic field traced by the OH masers ranges from 1.8 to 
16.3 mG, with a possible reversal. No OH masers were found associated with even the 
most prominent proplyds within 10 arcsec of 9 l Ori C. 
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1 INTRODUCTION 

The Orion A molecular cloud (L1641) is the nearest gi- 
ant molecular cloud (GMC) containing the nearest re- 
gions of massive star-formation. Because of its proximity 
it suffers little galactic extinction, and has been studied 
in far greater detail than other comparable regions in the 
Galaxy. The most recent formation of high-mass stars in the 
Orion GMC is within the Becklin-Neugebauer-Kleinmann- 
Low (BN/KL) region, about 1 arcmin Northwest from the 
Trapezium star 6 1 Ori C. The explosive nature of this event 
is demonstrated dramatically by the infrared H2 and [Fell] 
images of Allen & Burton (1993) and more recently Kaifu 
et al. (2000). 'Bullets' of matter are being ejected at hun- 
dreds of km s" 1 to form a bipolar cone of molecular 'fingers' 
whose axis is perpendicular to the hot core traced by NH3 
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emission (Wilson et al. 2000). These 'fingers' are tipped by 
Herbig-Haro (HH) objects whose shocked gas is visible at 
optical wavelengths (see Graham, Meaburn & Redman 2003 
for a recent association of these phenomena). 

At radio and millimetre wavelengths there are at least 
two molecular outflows, one associated with the optical and 
near-infrared features just described, and a second low- 
velocity outflow first detected through proper motions of 
H2O masers (Genzel et al. 1981). The powerful infrared 
source IRc2 was originally thought to be the source of the 
outflow, but more recent observations, in particular of SiO 
masers (e.g. Greenhill et al. 1998), suggest that radio con- 
tinuum source I of Menten & Reid (1995) is the more likely 
outflow source. Source I is significantly offset (by ~0'.'5) to 
the South of IRc2. 

Source I and the radio counterpart to BN have proper 
motions away from each other, corresponding to transverse 
speeds of 12 and 27 km s _1 respectively (Plambeck et al. 
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1995; Rodriguez et al. 2005). Here and elsewhere we as- 
sume a distance of 450 pc. The radio proper motions indi- 
cate that BN and source I were within 225 au of each other 
525±30 years ago. Tan (2004) has suggested that the large 
proper motion of BN is due to its ejection from the 1 Ori C 
system ~4000 years ago. Bally & Zinnecker (2005) on the 
other hand propose that the motions of BN and I are the 
result of a stellar merger ^500 years ago, which ejected BN 
and produced the energy powering the molecular outflow. 
Further analysis of the radio proper motion data by Gomez 
et al. (in press) shows that a third source n also has proper 
motions that project back to the position on the sky where 
BN and source I were close, ~500 years ago. This raises the 
possibility that n is the third star dynamically neccessary to 
eject BN from the multiple system. Distinguishing between 
these different possibilities is crucial to understanding the 
history of recent star-formation in the region. 

We have observed hydroxyl (OH) masers in the Orion- 
BN/KL region as part of an ongoing study of masers as- 
sociated with molecular outflows from massive young stars 
(Hutawarkorn & Cohen 2005, and references therein). OH 
masers have the potential to trace physical conditions, in- 
cluding magnetic fields which are detected and measured 
through Zeeman splitting and OH maser polarization. Pre- 
vious MERLIN OH observations by Norris (1984) detected 
80 1665-MHz masers distributed over a region of 10 arcsec, 
which were interpreted in terms of a rotating torus surround- 
ing IRc2. Johnston, Migenes & Norris (1989) detected 175 
masers spread over a region of 20 arcsec, with the majority 
in an E-W structure 14 arcsec in extent. The masers were 
found in clusters or 'clumps' that correlated with NH3 emis- 
sion. The present MERLIN observations more than double 
the number of OH masers known and show that the distri- 
bution is even more extensive in position and velocity (Sec- 
tion 3.1). The association of the masers with radio contin- 
uum sources is described in Sections 3.1 and 4.2, while the 
association with near- and mid-infrared sources is described 
in Section 3.3. 

Secondary targets in the same MERLIN observations 
were the 'proplyds' (O'Dell, Wen & Hu 1993) in the vicinity 
of 6 1 Ori C. These are the solar system-sized circumstel- 
lar envelopes discovered by Laques & Vidal (1979) which 
surround low-mass YSOs (Churchwell et al 1987; Meaburn 
1988; McCaughrean & Stauffer 1994; O'Dell & Wong 1996; 
Bally et al. 1998a). These stellar cocoons have dusty molec- 
ular disks with dense (> 10 6 cm~ 3 ) surfaces ionized by the 
intense flux of Lyman radiation from 6 1 Ori C (see Graham 
ct al. 2002). In these circumstances maser emission could be 
anticipated from the molecular disks which, if detectable, 
could be an important tool for investigating proto-planetary 
environments. Upper limits for such emission are reported 
in Section 3.5. 



2 OBSERVATIONS 

The OH observations were performed on 1998 April 27 and 
28 using seven telescopes of MERLIN (the Multi Element 
Radio Linked Interferometer Network): the 76-m Lovell 
Telescope and the Mk2 telescope at Jodrell Bank, and out- 
station telescopes at Pickmere, Darnhall, Knockin, Defford 
and Cambridge. The longest baseline was 218 km, giving a 



minimum fringe spacing of 0.16 arcsec. Observations covered 
all four ground-state OH transitions, namely 1612, 1665, 
1667 and 1720 MHz. A spectral bandwidth of 500 kHz was 
used, corresponding to a velocity range of 90 km s _1 . The 
spectral band was divided into 256 frequency channels, giv- 
ing a velocity spacing of 0.35 km s~ . The radial velocity at 
the centre of the 500-kHz band was +28 km s _1 with respect 
to the local standard of rest (lsr). Left-left and right-right 
circularly polarized signals from each pair of telescopes were 
simultaneously correlated to give LL and RR. 



The OH lines were observed in pairs, 1612 and 
1667 MHz for 6.5 hours on 27 April 1998 and 1665 and 
1720 MHz for 8 hours on 28 April 1998. Observations con- 
sisted of 5-min scans tracking the field centre RA(B1950) 
05 h 32 m 49 , :043, Dec. (B1950) -05°25'16'.'003 at each OH line 
frequency interleaved with 4-min scans on the nearby un- 
resolved phase calibrator source IC0539-057. Because of 
the faintness of the phase calibrator source it was neces- 
sary to observe it in 16-MHz bandwidth. Short ~ 30-min 
scans of the point source amplitude and bandpass calibra- 
tor 2134+004 were also made. 2134+004 was also used to 
provide the flux scale (based on comparisons with 3C286). 

The data processing and analysis procedures were car- 
ried out in B1950 coordinates, as described by Gasiprong, 
Cohen and Hutawarakorn (2002). Data were first edited, 
calibrated and corrected for gain-elevation effects using the 
Jodrell Bank d-programmes, and then passed into the AIPS 
software package. Within the AIPS package the data were 
further calibrated for all remaining instrumental and at- 
mospheric effects including the instrumental polarization. 
Self-calibration imaging of a bright pointlike reference chan- 
nel was employed to derive the final gain solutions for the 
1665-MHz data. These gain solutions were then applied to 
the other channels. The 1612- and 1667-MHz data sets had 
poorer signal-to- noise and could not be self-calibrated, while 
at 1720 MHz no emission was detected at all. Finally, maps 
in each circular polarization (RHC and LHC) were produced 
using CLEAN algorithms in AIPS, with a restoring beam 
of 150 mas x 150 mas. For each data cube the pixel size 
was 45x45 mas 2 and the total area mapped was 1024x1024 
pixels 2 (46x46 arcsec 2 ), centred on the position of IRc2. A 
second field of 512x512 pixels (23x23 arcsec 2 ) centred on 
the position of G 1 Ori C was also searched for possible emis- 
sion associated with the proplyds. 



The rms-noise in the final maps was 30-50 mjy beam , 
so only masers brighter than ~0.1 Jy could be detected. 
The masers were all unresolved at 150-mas resolution. Po- 
sitions of maser features were determined by fitting two- 
dimensional Gaussian components to the brightest peaks in 
each channel map and taking flux-weighted means across 
the channels showing emission from each particular feature. 
Final B1950 coordinates of masers were then converted to 
J2000 using the software package coco. The errors in the 
relative positions are typically 10 mas, while the absolute 
positions have a possible systematic error of 20 mas due to 
phase-transfer from the reference source. 
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Figure 1. Positions of OH masers (blue symbols, this paper), H2O masers (red symbols, Gaume et al. 1998), and compact continuum 
sources (green circles, Zapata et al. 2004) in Orion-BN/KL. The central region around IRc2 and radio source I is shown on an expanded 
scale in Fig. 2. The radio continuum sources BN, I and n are highlighted as filled green circles, and their proper motions are indicated 
schematically (data from Rodriguez et al. 2005 and Gomez et al. 2005). The direction to 8 1 Ori C is indicated by the dashed arrow. 
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Figure 2. Positions of OH masers (blue symbols, this paper), 
H2O masers (red symbols, Gaume et al. 1998) and SiO maser 
clusters (yellow symbols, Doleman et al. 1999) in the central re- 
gion of Orion-KL, near radio source I (Menten & Reid (1995). 
The position for source I corresponds to epoch 2000 (Rodriguez 
et al. 2000). The dashed ellipse corresponds to a ring with rota- 
tion axis at position angle —35° and inclination angle 58°to the 
line-of-sight, centred on radio source I (see Section 4.1). 



3 RESULTS 

3.1 Maser distribution in Orion-BN/KL 

We detected 11 masers at 1612 MHz, 430 masers at 
1665 MHz, 3 masers at 1667 MHz (LHC only) and no masers 
at 1720 MHz. The positions and velocities are given in Ta- 
ble 1. To facilitate comparison with other data, we give both 
the B1950 coordinates and the J2000 coordinates. The maser 
distribution is shown in Fig. 1 in J2000 coordinates. The OH 
masers are spread over a region 30 arcsec in extent, corre- 
sponding to ~ 13500 au. The radial velocities range from 
-13 to +42 km s _1 . The distribution is more extensive in 
both angular scale and velocity range than found in pre- 
vious work. However most of the emission comes from the 
dominant E-W structure noted by Norris (1984) and John- 
ston et al. (1989). 

One striking aspect of the new maps is that many of 
the OH masers lie in semi-continuous streams or arcs, of 
which the most prominent are labelled Stream A and Arc B 
in Fig. 1. Stream A stretches from RA (2000) 05 h 35 m 14:2, 
Dec. (2000) -05°22'25"to 05 h 35 m 14!5, -05°22'28", and in- 
cludes the maser 'clump' NW1 identified by Johnston et 
al. (1989). The radial velocities are approximately constant, 



around +21±2 km s" 1 . Arc B stretches from 05 h 35 m 14T, 
-05°22'28"to 05 h 35 m 14=3, -05°22'28", and includes the 
maser clumps NW2 and NW3 of Johnston et al. (1989). 
The radial velocities cover a wide range, from ~-7 km s _1 
in the East to ~+10 km s _1 in the West. These streams or 
arcs are large-scale features, with angular sizes of ~5 arc- 
sec that correspond to ~2250 au. We argue that they trace 
large-scale shocks or ionization fronts (Section 4.2). 

The distribution of H2O masers from Gaume et al. 
(1998) is also shown in Fig. 1. The OH and H2O masers 
have complementary distributions, with almost no overlap. 
An expanded plot of the crowded central region is shown in 
Fig. 2, revealing the same situation on the small scale. The 
H2O masers in this region are the so-called 'shell' masers 
that were not detected in the VLBI measurements by Genzel 
et al. (1981). There is an elliptical zone of avoidance centred 
on source I, within which there are very few OH masers, but 
the H2O 'shell' masers and SiO masers are found. The H2O 
'shell' masers lie inside the dashed ellipse in the figure, while 
nearly all of the OH masers lies outside the ellipse. The po- 
sition angle and inclination angle of the ellipse were chosen 
to match the model fitted to the OH masers in Section 4.1. 

Fig. 2 also shows the positions of the SiO maser clusters 
mapped by Doeleman, Lonsdale & Pelkey (1999). The SiO 
masers cluster tightly around the position of the radio source 
I (Menten & Reid 1995), with the H2O 'shell' masers sur- 
rounding them. The relative location of SiO, H2O shell and 
OH masers at increasing distances from the central source I 
is similar to that found for these masers in the circumstellar 
envelopes of post-AGB stars, which has a natural expla- 
nation in terms of excitation by a powerful central source 
(Chapman & Cohen 1986). 

The positions of compact radio continuum sources in 
the region are plotted in Fig. 1 as open green circles (data 
from Zapata et al. 2004), with sources I, n and BN high- 
lighted as filled green circles. Apart from source I, which is 
at the centre of the SiO and H2O 'shell' maser distributions 
and the OH zone of avoidance, there is another association 
of OH masers with a compact radio continuum source that is 
noteworthy, namely a striking association between Stream 
A and BN. The stream points directly towards BN (away 
from IRc2), in the direction of BN's proper motion (shown 
schematically by the arrow in Fig. 1). This is discussed fur- 
ther in Section. 4.2. Finally we note a possible weaker associ- 
ation of source D of Menten & Reid (1995) (source 21 of Zap- 
ata et al.), near RA (J2000) = 05 h 35 m 14=90, Dec. (J2000) = 
-05°22'25'.'4, with a cluster of OH masers, corresponding to 
'clump NET of Johnston et al. (1989). The radial velocities 
are centred near +8 km s _1 . 

Further study of Fig. 1 reveals that the OH 1612-MHz 
masers are distributed differently from the mainline 1665- 
and 1667-MHz masers. Most of the OH 1612-MHz masers 
are found a long way from the centre, in the outskirts of 
the OH maser region, and they are spatially separated from 
the other OH masers. Their far-flung distribution is strik- 
ingly different from that of the other OH masers, but similar 
to that of the H2O masers. This similarity is reinforced by 
kinematic similarities. 
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Figure 4. Position-velocity plot showing declination (J2000) and radial velocity for the OH masers (blue symbols, this paper) and the 
H2O masers (red symbols, Gaume et al. 1998) in Orion-BN/KL. 
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3.2 Maser kinematics 

The complex kinematics of the Orion-BN/KL masers are 
summarized in the position-velocity diagrams Figs. 3 and 4. 
Fig. 3 shows RA vs. VJ sr for the OH and H2O masers. There 
is a clear contrast between the two species: the OH masers 
show a dominant rotational pattern, in which the radial ve- 
locities drop from +20 km s _1 in the East to km s _1 in 
the West, while the H2O masers show a dominant expan- 
sion pattern, with highly red-shifted or blue-shifted masers 
equally likely to occur to the East or West. This is con- 
sistent with the well known expansional proper motions of 
these masers (Genzel et al. 1981; Gaume et al. 1998). The 
OH rotational pattern is clearest at right ascensions greater 
than05 h 35 m 14!5, and more disturbed at right ascensions less 
than this. Stream A has a velocity greatly different from the 
rotational pattern, as indicated in Fig. 3. We note that the 
radial velocity of Stream A, +21 km s _1 , is the same as that 
of BN and its circumstellar nebulae within 25 au (Scoville 
et al. 1983). 

Fig. 4 shows a Dec-Vi sr plot. The OH masers trace a 
dominant inclined elliptical pattern, centred on Dec. (2000) 
-05°22'29"and V lsl = +6 km s" 1 . Radial velocities rise 
abruptly from 0-10 km s _1 in the South to +20 km s _1 
to the North of Dec. (2000) -05°22'29". The 'hole' in the 
centre of the Dec-Vi sr ellipse suggests that in addition to 
rotation there is an expansional component of motion of 
similar magnitude. The locations of OH and H2O masers in 
this diagram are again complementary. 

The OH 1612-MHz masers are kinematically and spa- 
tially distinct from the 1665- and 1667-MHz masers. This is 
particularly clear in Fig. 4, where they are well separated 
in both declination and velocity, appearing at the extremes 
of both coordinates, well separated from the rotating and 
expanding ring signature of the mainline masers. The kine- 
matic pattern is similar to but even more extensive than 
that of the H2O masers, suggesting that outflow rather than 
rotation is the dominant motion. 

The OH distribution and kinematics are modelled in 
Section 4.1. 

3.3 Comparison with Infrared Data 

The distribution and kinematics of the dominant OH emis- 
sion show a degree of symmetry about the position of IRc2 
and radio source I, and a radial velocity of ~8 km s _1 , as 
noted by previous authors cited in the introduction. The 
OH masers encompass IRc2, 3, 6 and 7. The availability of 
high resolution infrared data allows us to examine in more 
detail the correspondences of the OH masers with the pow- 
erful mid-infrared sources and with near-infrared bullets in 
the molecular outflow. 

3.3.1 Comparison with 12.5-fim Keck and 11.7-fxm 
Gemini Data 

Shuping, Morris & Bally (2004) have published a high reso- 
lution (0'.'38) mid- infrared image of the Orion BN/KL region 
obtained with the Keck I telescope at 12.5-/im. Fig. 5 is a 
composite plot showing the OH masers (this paper) superim- 
posed on Fig. 5 of Shuping et al. (2004). This shows several 
surprising results. Firstly, OH masers cluster around IRc2, 



but they avoid the actual source position (relative coordi- 
nates (+5", -7") in this diagram), which lies in the zone of 
avoidance noted in Section 3.1. Secondly, Stream A points 
directly away from IRc2 towards the dominant mid-infrared 
source BN (at the origin), at position angle ~45°. Stream A 
runs through a region devoid of mid-infrared emission. The 
direction from IRc2 to BN is parallel to the proper motions 
recently measured for BN and radio source I by Rodriguez 
et al. (2005). It is also the projected direction of the large- 
scale molecular outflow. Thirdly, Arc B curves parabolically 
around the extended source IRc6 (near (+l",-4") in this 
figure), in a manner that suggests a close physical associa- 
tion. This confirms and extends the result of Johnston et al. 
(1989). 

Smith et al. (2005) have published a high-resolution 
11.7-/im mosaic image of the inner Orion nebula obtained 
with the Gemini South telescope. The only close correlation 
(better than 2 arcsec) of any of the OH 1.6-GHz masers 
listed in Table 1 with the 11.7-/im point sources listed by 
Smith et al. (2005) is the 1612L maser number 10, at RA 
(J2000) = 05 h 35 m 14:6533, Dec. (J2000) = -05°22'507022, 
which has a radial velocity of +18.9 km s _1 . The corre- 
sponding IR source is at 05 h 35 m 14:67, -05°22'49'.'5. Smith 
et al. (2005) describe the latter as 'an embedded IR source 
with no optical ID'. None of the 27 known proplyds that 
correlate with 11.7-/u,m point sources in Smith et al. (2005) 
are identified here as maser sources. 

3.3.2 Comparison with 2.12-fxm Subaru Data 

The locations of the OH masers relative to the well-studied 
molecular outflow are shown in Fig. 6. The infrared H2 2.12- 
fim images presented here were taken with the Subaru 
Telescope as part of its astronomical first light programme 
(Kaifu et al. 2000). The data were taken on 1999 January 
11, 13 and March 6 using the Cooled Infrared Spectrograph 
and Camera for the OH-suppression spectrograph, CISCO 
(Motohara et al. 1998). 

The 2.12-/im image shown in Fig. 6 is a 150" x 
150"section of the original 5' x 5' mosaic (provided by Masa 
Hayashi). In order to overlay the OH maser positions onto 
the Subaru image, it was first necessary to astrometrically 
register the image with the J2000 coordinate system. This 
was done using > 600 stars from the catalogue of Hillen- 
brand &: Carpenter (2000) and the Starlink GAIA software 
package. A computer programme was then written to read 
the B1950 OH maser positions and convert them to J2000 
coordinates using the Starlink Astrometry Library and pro- 
duce a catalogue that was readable by the GAIA software 
package. The positions were then overlaid using GAIA. 

The 1665- and 1667-MHz masers (shown by dot and 
cross symbols respectively) are concentrated towards the 
highly obscured IRc2 region. In general they show little cor- 
respondence with features in the near-IR map. However the 
1612-MHz masers (+) correspond with fingers of H2 emis- 
sion in the molecular outflow. This association is seen more 
clearly in Fig. 7, which shows the central region on an ex- 
panded scale. The 1612-MHz masers show a strong associ- 
ation with fingers of H2 emission in the molecular outflow, 
particularly in the South, near RA (J2000) 05 h 35 m 14:7, Dec. 
(J2000) -05°22'46", where 6 (half) of the 1612-MHz masers 
are found. Three 1665R masers (numbers 63, 64 and 76 in 
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Figure 5. Positions of OH masers (black symbols as in Fig. 1) overlaid on the 12.5-/im image (greyscale) taken with the Keck Telescope 
by Shuping et al. (2004), adapted from their Figure 5. The open circles show the positions of OH maser clumps identified by Johnston et 
al. (1989). Red and blue dots indicate H2O masers from Gaume et al. (1998) and Genzel et al. (1981), with red indicating radial velocity 
greater than +5 km s — 1 and blue indicating radial velocity less than +5 km s . The green contours show NH3 (J,K) = (4,4) emission 
associated with the 'hot core', from Wilson et al. (2000). Offsets are relative to the position of the BN object: RA (J2000) = 05 h 35 m 14?117, 
Dec. (J2000) = -05°22'22'/9. IRc2 is near (+5", -7") and IRc6 is near (+l",-4"). 



Table 1) are found nearby, and may also be associated. This 
is the first time that clear counterparts to interstellar OH 
masers have been seen. It is not clear, however, why these 
particular H2 emission fingers should have OH counterparts, 
while many others do not. 

3.4 Zeeman Splitting and Magnetic Fields 

The data were searched for possible Zeeman pairs of opposite 
circular polarization, with radial velocities differing by more 
than 0.8 km s" 1 (to exclude linearly polarized features) and 
lying within 30 mas of each other. The results are given in 
Table 2. We estimate the chance of a false association to be 
no more than 4 percent. 

The magnetic fields range from -3.5 mG (directed to- 
wards us) to +16.3 mG. In Fig. 8 we plot the magnetic field 



values as a function of position on the sky, together with 
values from previous interferometric OH Zeeman measure- 
ments (Hansen & Johnston 1983; Norris 1984; Johnston et 
al. 1989). The present data suggest that in addition to the 
overall field of 1-3 mG directed away from us, there may 
also be regions of higher field strength (Z4), and that the 
field direction may reverse (Z3). The field reversal is seen 
at only a single point, and not on the larger-scale that is 
found in other bipolar outflow sources (Blaskiewicz et al. 
2005; Hutawarakorn & Cohen 2005, and references therein). 
More data are needed to confirm these findings. 

3.5 Search of the Proplyd Region 

A search of a 23 arcsecx23 arcsec region around 6 1 Ori C was 
made, but no OH masers with peak flux densities >0.1 Jy 
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Figure 6. Positions of the OH masers in Orion-BN/KL, overlaid on the infrared 2.12-/im H2 image taken with the Subaru Telescope 
by Kaifu et al. (2000), showing the relationship of the masers to the outflow (Section 3.3). The 1665-MHz masers (.) and the 1667-MHz 
masers (x) lie in an expanding and rotating disc or torus centred on IRc2 and source I (Section 4.1), while the 1612-MHz masers (+) 
have a more widespread distribution. The area in the box is shown on an expanded scale in Fig. 7. 



were found, even near the most prominent proplyds (LV 1- 
6). The densities in these regions have been estimated to be 
10 5 -10 6 cm (Section 4.4), which are typical of OH maser 
sources, while the temperatures are estimated to be 200- 
350 K (Hayward, Houck & Miles 1995), which are somewhat 
higher than usually encountered in OH maser sources (e.g. 
Cragg, Sobolev & Godfrey 2002). The non-detection of OH 
masers might be explained by a low OH column density, 
which could arise through chemical evolution (e.g. Rodgers 
& Charnely 2001). 



4 DISCUSSION 

4.1 Kinematic Model for OH 1665-MHz Masers 

For our kinematic modelling we used only the 1665-MHz 
masers. The celestial coordinates and radial velocities of 
the OH masers give incomplete information on the three- 
dimensional source geometry and motions. The OH maser 
positions are essentially the x- and y-components of each 
maser spot, but the z-coordinate is unknown. The radial 
velocity of a maser is related to the unknown velocity field 
and the (x,y,z) coordinates of the maser. We modelled the 
source using a simple kinematic model which was fitted to 
the data using the Control Random Search technique (Price 
1976). The model enables us to find the z-coordinate for 
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Table 2. Possible OH Zeeman pairs in Orion. 



Zccman 


Transition 


Vi sr 


RA (1950) 


Dec. (1950) 


RA (2000) 


Dec. (2000) 


Spk 


B 


Pair 




km s 


05 h 32 m s 


-5°:24': " 


05 h 35 m s 


-5°:22': " 


(Jy) 


mG 


7,1 


1665L 


21.98 


4fi 8077 


19.522 


1 4 9^0 




0.18 


4-1 8 




1665R 


23.06 


46.8102 


19.504 


14.2855 


26.047 


0.10 




z 2 


1665L 


18.13 


47.3239 


21.992 


14.7984 


28.572 


0.10 


+3.6 




1665R 


20.24 


47.3210 


22.020 


14.7954 


28.600 


0.38 




z 3 


1665L 


15.67 


46.7458 


19.896 


14.2210 


26.434 


0.14 


-3.5 




1665R 


13.56 


46.7487 


19.915 


14.2239 


26.454 


0.22 




z 4 


1665L 


11.59 


46.9193 


23.131 


14.3934 


29.682 


0.26 


+16.3 




1665R 


1.96 


46.9203 


23.109 


14.3944 


29.660 


0.10 




z 5 


1665L 


0.30 


46.8021 


23.549 


14.2761 


30.092 


0.58 


+3.6 




1665R 


2.43 


46.7989 


23.565 


14.2729 


30.108 


1.48 






15 14 

Right Ascension (J2Q00) 



5:35:13 



Figure 7. Expanded view of the central region of Fig. 6, showing 
the positions of OH masers overlaid on the infrared 2.12-/im H2 
image taken with the Subaru Telescope by Kaifu ct al. (2000). The 
1612-MHz masers (+) arc associated with fingers of H2 emission 
in the molecular outflow. 
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Figure 8. Magnetic field measurements in the Orion-BN/KL re- 
gion from OH maser Zeeman splitting (present paper, plus ref- 
erences indicated in the legend). The locations and schematic 
proper motions of the continuum sources BN, I and n are shown 
for reference, as in Fig. 1. 



each maser from the radial velocity and the velocity field, 
and so produce a three-dimensional view of the OH maser 
distribution (Fig. 9). 

We assumed, for modelling purposes, that IRc2 is at 
the kinematic centre. However the results would not be af- 
fected if radio source I were taken as the centre. We es- 
timate that the kinematic centre of the OH masers is un- 
certain to ~1 arcsec. The masers were assumed to have a 
uniform expansion velocity V cxp away from IRc2, centred 
on an unknown radial velocity Ko, plus solid body rotation 
corresponding to Kot at 1-arcsec separation from IRc2. The 
rotation axis is tilted by an angle a about the y-axis and 
an angle (5 about the x-axis from the x — z plane (follow- 
ing Reid et al. 1988). The parameters Kx P , Ko, Vrot, a. and 



f3 were the five unknowns. For each maser and each posi- 
tion in the five-dimensional parameter space, we compared 
the observed and calculated radial velocity V\ BI and found 
the z-position where the absolute value of the residual was 
minimized (allowing the maser to lie anywhere along the z- 
axis, with a distance corresponding to 30 arcsec of IRc2). 
The five-dimensional parameter space was searched, within 
a reasonable domain, using a random number generator to 
sample the parameters. The position in five-parameter space 
where the sum of the (absolute) velocity residuals is mini- 
mized was thereby located. The best fitting model has the 
parameter values given in Table 3. The angles correspond 
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Figure 9. Three-dimensional view of the 1665-MHz OH masers in Orion-KL, according to the kinematic model that was fitted (Sec- 
tion 4.1). The x— and y— axes correspond to RA and Dec, while the z— axis is in the line-of-sight direction. Box A shows the view in the 
x — y (RA-Dec) plane looking in the +z— direction, box B shows the view looking in the +x— direction (increasing RA), box C shows the 
view looking along the rotational axis of the best-fitting model, and box D shows the view looking down the y— axis (decreasing Dec). 



to a rotation axis inclined at 58°to the line-of-sight with a 
position angle on the sky of -35° . A disc at this orientation 
is plotted in Fig. 2, where it fits neatly inside the 'hole' in 
the OH maser distribution that is centred on source I. 

Using the kinematic model we constructed views of 
the three-dimensional distribution of OH 1665-MHz masers, 
which are shown in Fig. 9. The masers lie in an irregu- 
larly filled torus, at radial distances ranging from 430 au 
to 13200 au, with a mean radius of 6200 au. There is a well 
defined inner cavity, of radius ~1300 au. This cavity cor- 
responds spatially with the SiO flared disc region mapped 
by Wright et al. (1995, 1996), which has a radial velocity 
range of —10 to +20 km s _1 that is roughly consistent with 
the range of —12 to +30 km s _1 given by our kinematic 



Table 3. Parameters of kinematic model. 



Parameter 


Search Range 


Best fit 


Vcxp 


15-32 km s" 1 


21.0+3.5 km s" 1 




8-10 km s" 1 


9.0±0.5 km s" 1 


Vrot 


1—5 km s — 1 


2.9±0.9 km s" 1 


a 


0-180° 


46°±22° 


P 


0-90° 


48°±19° 



model. The distribution in radial distance R" from the ro- 
tation axis is plotted in Fig. 10. The distribution in the 
^''-direction, parallel to the rotation axis, has a full width 
to half maximum of 6000 au, and a total extent of 12000 au. 
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Figure 10. Radial distribution of the 1665-MHz OH masers 
about the disc rotation axis, determined using the kinematic 
model (Section 4.1). 



The best-fitting expansion velocity is similar to that of 
the 'low- velocity outflow' seen in H2O masers (Genzel et al., 
1981). Solid body rotation becomes the dominant motion be- 
yond a distance of 3240 au from the centre (a distance corre- 
sponding to 7.2 arcsec). Rotation is therefore the dominant 
motion for most of the system of OH masers, with how- 
ever a significant component of expansion. The rotational 
period in our model is 760 yr, while the expansion timescale 
varies from ~300 yr for the inner edge of the maser cavity 
to ~3000 yr for the masers most distant from the expansion 
centre. 

Most of the Orion nebula HH-objects are found at the 
tips of the H2 'fingers' in the bipolar outflow shown in Fig. 6 
and appear to have their origin in a newly created massive 
star that gives rise to the infrared source IRc 2 in the 10 5 
L© Orion-BN/KL nebula in the OMC-1 cloud core. It has 
also been suggested that a minority of the prominent HH- 
objects have their origin in the neighbouring star forming re- 
gion OMC-1S (Smith et al. 2004). However, excluding these 
few, Doi, O'Dell & Hartigan (2002) have used HST proper- 
motion measurements of the global expansion of the system 
of Orion HH-objects to demonstrate that the bipolar out- 
flow is of Hubble type, i.e. the flow velocity is proportional to 
distance from the expansion centre. The bipolar flow there- 
fore consists of 'bullets' all ejected around a 1000 yr ago 
in a single explosive event. Thus the expansion timescales 
determined for the OH masers bracket this H2 timescale. 
These results suggest that the OH maser disc acquired its 
expansional motions in the same event that produced the 
H2 bullets and the HH-objects found at the tips of the H2 
'fingers' in the bipolar outflow shown in Fig. 7. 

The orientation of the rotation axis determined here 
agrees with that of the molecular outflow traced in CO and 
H 2 (e.g. Chernin & Wright 1996; Kaifu et al. 2000). It also 
coincides with the direction of the proper motions of source 
I and BN, and is close to the large-scale magnetic field di- 



rection inferred from far-infrared and submillimetre polar- 
ization (Schleuning 1998 and references therein). 

We note that the velocities of order 20 km s _1 im- 
plied by our kinematic model correspond to proper motions 
^9 mas yr -1 , which could be measured by MERLIN within 
ten years. 



4.2 The Nature of Stream A and Arc B 

The major streams or arcs found in Orion-BN/KL have an- 
gular sizes greater than 5 arcsec, corresponding to projected 
sizes of at least 2250 au, which are almost unprecedented 
for interstellar OH masers. Stream A has no obvious coun- 
terpart in the infrared. However it points in the direction 
of the dominant infrared source BN, in the direction of the 
radio proper motions of BN and source I. Stream A has the 
appearance of a vapour trail that has formed in the wake of 
the runaway star BN. From the proper motion of BN and 
the projected separation from the Northern tip of Stream 
A we can infer a timescale of 200 yr. This may correspond 
to a cooling time, a formation time for OH, or a time for 
population inversion to be established following the passage 
of BN. 

Fig. 1 shows that the proper motions of BN and sources 
I and n project back to a position on the sky near the base 
of Stream A, in a region largely clear of OH masers. This 
position is significantly displaced ~4 arcsec Northwest of the 
centre of the OH maser torus. It has been suggested that BN 
and sources I and n originally belonged to a multiple stellar 
system that disintegrated ~500 yr ago. The disintegration 
could have due to a close dynamical interaction, as suggested 
by Gomez et al. (in press), or it could have been due to a 
merger event that also produced the molecular outflow, as 
suggested by Bally & Zinnecker (2005). Alternatively, Tan 
(2004) has proposed that BN was ejected from the 9 Ori 
C system ~4000 yr ago. Our data do not distinguish un- 
ambiguously between these possibilities. However the fact 
that the large-scale OH maser torus is centred on source 
I slightly favours Tan's scenario, with BN being a passing 
runaway from the 8 1 Ori C system. 

The other large-scale feature Arc B encloses IRc6, which 
is one of the cooler mid-IR sources, with polarization that is 
consistent with external illumination and with no evidence 
of embedded bright stars (Shuping et al. 2004 and references 
therein). This suggests that Arc B might be associated with 
a photodissociation zone around IRc6. In further support of 
this there is strong emission from CH3CH at the position 
of IRc6, implying densities of at least 10 6 cm -3 (Wilner, 
Wright & Plambeck 1996). 

The large-scale maser structures found in Orion-BN/KL 
are similar in scale to OH maser structures recently re- 
ported in W3(OH). Wright et al. (2004a,b) noted several 
arcs of 18-cm masers, which they interpreted as propagating 
shocks. Harvey-Smith & Cohen (2005) found a low bright- 
ness filament of excited OH 4765-MHz emission stretching 
for ~2200 au and clearly related in its morphology to some 
of the ground-state OH arcs at 18-cm wavelength. The struc- 
tures seen in Orion are similar in linear scale. The individual 
maser spots in Stream A and Arc B are well separated on 
the sky for the most part, but are nevertheless likely to be 
simply the high-gain cores of an extended stream of maser 
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emission. The typical flux densities of ~0.1 Jy correspond 
to brightness temperatures of at least ~3xl0 K. 

All these large-scale features in Orion and W3(OH) are 
relatively weak and therefore below the sensitivies of earlier 
surveys, for the most part. Deeper observations of other well- 
known OH maser sources are needed to establish just how 
common such features are. 



4.3 The 1612-MHz Masers 

The OH 1612-MHz masers have a more widespread distri- 
bution than the other OH masers, and different kinematics 
(Section 3.2). This suggests that they trace a region with 
different physical conditions. Gray, Field & Doel (1992) give 
the following conditions for strong 1612-MHz masers: molec- 
ular hydrogen density riH 2 =6 x 10 6 cm -3 , «oh 60 cm -3 , gas 
kinetic temperature 21=30 K, dust temperature at maser 
site Td 30 K, external radiation field T x 80 K and veloc- 
ity shift AV^ 2.0 kms -1 . This is cooler than typical for OH 
masers, but is consistent with the location of these partic- 
ular 1612-MHz masers further than normal from the main 
source of infrared luminosity. The velocity gradient needed 
for strong maser action is also consistent with the association 
of these masers with shocked gas in the molecular outflow 
(Section 3.3.2). 

We note a further possible positional association be- 
tween the 1612-MHz masers and the methanol 6.7-GHz 
masers recently discovered by Voronkov et al. (2005). In 
Fig. 11 we show the positions and velocities of the two 
species. There is a clear region of overlap centred around 
RA (J2000) = 05 h 35 m 14:5, Dec. (J2000) = -05°22'45'.'9 and 
Visr= +7.5 km s _1 . The 6.7-GHz positions have errors of 
~2 arcsec. It will be important to make 6.7-GHz measure- 
ments of higher precision to examine this correspondence 
more closely. The positions of 25-GHz masers from John- 
ston et al. (1997, 1992) are also shown in the figure for com- 
pleteness. Only one of these falls in the region of interest. 
We note that there is no association with the OH mainline 
masers. This is consistent with the 25-GHz masers being 
class I masers, which are thought to be collisionally pumped 
(cf. Cragg et al. 2002). 

The 6.7-GHz maser transition is the prototype class II 
maser, usually thought to be radiatively pumped. The 6.7- 
GHz masers in Orion are of therefore of interest because of 
their apparent association with 25-GHz masers. Voronkov 
et al. have considered the pumping requirements for the co- 
existence of both types of methanol maser and find that 
both types can occur simultaneously at low temperature 
(~60 K) and low molecular hydrogen density (~10 5 cm -3 ). 



These conditions are not too dissimilar to those needed for 
strong 1612-MHz masers, as given earlier. In summary, it ap- 
pears that in Orion we may have the first examples of both 
methanol class II and OH ground state masers located far 
from the main source of infrared luminosity, and associated 
instead with the interaction between the molecular outflow 
and the surrounding gas. 



4.4 Proplyd region 

The proplyds originally designated LV 1-6 within the Orion 
Trapezium cluster are those most likely to exhibit maser 
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Figure 11. Positions and velocities of OH 1612-MHz masers (this 
paper) compared with those of the newly discovered methanol 6.7- 
GHz masers (Voronkov et al. 2005). Positions of 25-GHz methanol 
masers from Johnston et al. (1997) and Johnston et al. (1992) 
are also shown. The radial velocities at 25 GHz are confined to a 
narrow range of 7-10 km s — 1 . 



emission for they are irradiated by the intense field of the 
06.5 star 6 1 Ori C. Of these LV 1, now resolved as a binary 
proplyd (168-326 NW & SE) by O'Dell & Wen (1994), and 
LV 2 (proplyd 167-317 of O'Dell & Wen 1994) have been 
most intensively observed (Graham et al. 2002; Henney et 
al. 2002 respectively). 

The 8xl0 14 cm radius LV 2 was shown by Henney et 
al. (2002) to be irradiated by a flux of 3xl0 4 cm" 2 s" 1 
Lyman photons. These produce an ionised skin with a mea- 
sured electron density n e = 2.6 xlO 6 cm -3 on the surface of 
a neutral disk or cocoon which envelopes the low mass YSO 
(Meaburn 1988). From 1.3-mm interferometric observations 
Bally et al. (1998b) deduce an upper limit to the mass of 
molecular hydrogen in this neutral component of ~0.015 Mq 
to give an H2 density of n < 4.2 xlO 9 cm -3 . On the other 
hand, simple pressure balance between the ionised and neu- 
tral gas with temperatures of T e = 10 4 K and T = 10 3 K 
would give n w 5.2xl0 9 cm -3 . 

OH masers have yet to be found in regions of purely 
low-mass star-formation. However, it is this dense, circum- 
stellar molecular gas around proplyds, stimulated by pho- 
tons from 6 1 Ori C, that is the potential source of any maser 
emission. We would expect to find OH masers at a distance 
of ~7xl0 17 cm from an 06.5 star (cf. Baart & Cohen 1985). 
However, if Keplerian motion around the YSO occurs in this 
neutral cocoon it will have a rotational velocity of 4 km s _1 
at its outer edge with higher velocities towards the centre 
(proportional to radius -1 / 2 ). In these circumstances the ve- 
locity coherence required for strong maser amplification will 
occur over a path length that is much smaller than the over- 
all radius of the proplyd. This, together with possible chem- 
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ical effects (Section 3.5), could explain our non-detection of 
OH masers in this region. 



5 CONCLUSIONS 

The distribution of OH masers in the Orion-BN/KL region 
is far more extensive than previously realized, covering a re- 
gion of 30 arcsec extent and a radial velocity range from -13 
to +42 km s _1 . The bulk of the emission can be modelled in 
terms of a rotating and expanding torus, centred on IRc2 or 
radio source I, with an inner cavity of ~1300 au radius. The 
rotation axis has the same position angle and inclination to 
the line-of-sight as the molecular outflow and the large scale 
magnetic field inferred from mm- and submm-polarization 
(Section 4.1). The dynamical timescale is similar to that of 
the explosive event that produced the widespread shocked 
H2 emission. It is likely that the OH masers trace the inter- 
action between the low-velocity molecular outflow and the 
molecular hot core, and that they acquired their expansional 
motions in the same event that produced the outflow. 

Of particular interest is a string of masers, Stream A, 
at ~21 km s~\ that extends at position angle ~45°between 
IRc2 and BN, in the direction of the radio proper motions of 
these two dominant sources. We suggest that Stream A may 
have appeared, like a vapour trail, in the wake of the run- 
away star BN. The proper motions of BN and sources I and 
n project back to the base of Stream A (Fig. 1), ~4 arcsec 
Northwest of the centre of the OH maser torus, a position 
that is largely devoid of masers. 

The 1612-MHz masers have a more widespread distri- 
bution than the other OH masers, with kinematics that are 
more like those of the H2O masers associated with the out- 
flow. Many of these 1612-MHz masers are spatially associ- 
ated with fingers of shocked H2 emission (Fig. 7). These OH 
masers are thought to require relatively low gas and dust 
temperatures for their inversion (Section 4.3). Apart from 
the OH 1612-MHz masers, the other OH masers have com- 
plementary distributions and kinematics to the H2O masers 
(Figs. 1-4), with essentially no overlap. OH is also spatially 
and kinematically distinct from the class I methanol 25-GHz 
masers. A possible correspondence between OH 1612-MHz 
and class II methanol 6.7-GHz masers reported by Voronkov 
et al. (2005) requires further study at higher angular reso- 
lution. 

The magnetic field strength in the OH maser regions 
ranges from 1.8 to 16.3 mG, with a possible field reversal 
(Section 3.4). 

A search for OH masers associated with the proplyds 
around 9 1 Ori C yielded only upper limits, which may indi- 
cate a low OH abundance in this more evolved region. 
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OH maser disc and outflow in Orion-BN/KL 



Table 1. OH 1.6-GHz Masers in Orion-BN/KL. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1665-MHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 



Transition Feature 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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Table 1. OH 1.6-GHz Masers in Orion-BN/KL, continued. 
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